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The time-and-frequency resolved nonlinear light scattering (NLS) signals from a time evolving 
charge distribution of valence electrons prepared by impulsive X-ray pulses are calculated using a 
superoperator Green's function formalism. The signal consists of a coherent ~ iV 2 -scafing difference 
frequency generation and an incoherent fluorescence ~ iV-scaling component where N is the number 
of active molecules. The former is given by the classical Larmor formula based on the time-dependent 
charge density. The latter requires additional information about the electronic structure and may 
be recast in terms of transition amplitudes representing quantum matter pathways. 
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I. INTRODUCTION 

Newly developed attosecond X-ray sources can excite 
molecular electronic states impulsively, paving the way 
for novel spectroscopic probes of electronic structure and 
correlations [1-10]. Thanks to their broad bandwidth, 
such pulses can prepare molecules or molecular ions in 
non-stationary superpositions of electronic states [11]. 
These states can be probed by various nonlinear optical 
techniques, such as photoelectron spectroscopy [12, 13]. 

Traditional X-ray experiments [14], including X-ray 
absorption near-edge structure (XANES), resonant and 
nonresonant X-ray emission spectroscopy (RXES and 
NXES), Auger electron spectroscopy (AES), and X-ray 
diffraction [15] (XD), are mostly related to single-particle 
characteristics of the many-body ground-state wave func- 
tion. Quantum coherence (i.e., the phase in a superpo- 
sition of states), which is accessible by nonlinear spec- 
troscopy, does not play a role in these techniques. Co- 
herent X-ray sources can look into these quantum effects. 
In particular, the large bandwidth (10 eV for an 100 as 
pulse) can be used to create electronic coherences which 
provide a higher level of information about orbitals that 
goes beyond the charge density [16]. Nonlinear effects 
in the X-ray regime have long been observed in the fre- 
quency domain [17], including parametric downconver- 
sion [18-20] hard X-ray frequency doubling [21] and two- 
photon X-ray fluorescence [22]. Combined optical/X-ray 
nonlinear techniques can monitor excitations of optically 
excited states [18, 23-25]. All-X-ray nonlinearities such 
as pump-probe have been demonstrated as well [11, 19- 
21, 26-28]. 

Classically, spontaneous emission is related to the ac- 
celeration of charges [13]. In the semiclassical approach 
(classical field/quantum matter) [29] the emission spec- 
trum is thus calculated from a time dependent charge 
density that acts as a source. No other information about 
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FIG. 1. (Color online) (a) - setup for the X-ray induced NLS, 
(b) - level scheme of the model system, go,..., gi represent 
valence states, / are core excited states. 



the matter is needed. This results in the Larmor for- 
mula for the power emitted by a radiating charge (Eq. 
15) [13, 30]. In a more rigorous description where both 
field and matter are treated quantum mechanically [31] 
we must work in the joint field/matter space and cal- 
culate the signal perturbatively in the radiation-matter 
coupling. Radiation back-reaction is then included im- 
plicitly, through the joint dynamics of field and matter. 

In the present work we study the nonlinear light scat- 
tering (NLS) caused by the time evolving superposition 
of electronic excitations. In order to create such super- 
positions the coherent X-ray source with broad (~10 eV) 
bandwidth is required. We use the many-body Green's 
function formalism [32] to study NLS corresponding to 
transitions between valence excitations triggered by stim- 
ulated X-ray Raman process [10, 33, 34] (see Fig. 1). 
Note, that NLS results from the single pulse used for 
excitation; there is no probe pulse to simulate the emis- 
sion. The resulting time and/or frequency gated signals 
[35, 36] provides direct information on valence electron 
motions, coherences and correlations [10, 37]. We iden- 
tify two possible mechanisms for this process to quadratic 
order in X-ray pulse intensity: difference frequency gen- 
eration (DFG) and stimulated Raman induced fluores- 
cence (SRIF). The former is a coherent process associated 
with long range valence state coherence, which can be 
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FIG. 2. (Color online) (a) loop diagram for the bare incoher- 
ent signal in a gated measurement caused by a single molecule 
a, (b) coherent signal generated by a pair of molecules a ^ f3. 
Eqs. (2) - (3) (coherent) and Eqs. (8) - (9) (incoherent) 
can be read off these diagrams (for diagram rules see Ref. 
[38]). The shaded area represents an excitation by an arbi- 
trary sequence of pulses, which prepares the molecule in a 
superposition state. 



described in terms of the time-dependent charge density 
created by the Raman process i.e. the classical Larmor 
formula. It is highly directional (phase matched) and 
scales as N(N — 1) with the number of active molecules. 
The latter is an incoherent process coming from excited 
state populations which requires more detailed informa- 
tion about the excited states than the simple charge den- 
sity. It produces an isotropic signal and scales as N. 
For small samples, the coherent process is dominated 
by a short range coherence coming from the neighbor 
molecules that gives rise to an isotropic emission and 
scales linearly with N. As an example we compare both 
contributions for cysteine excited at the S, O and N core 
transitions. 

The paper is organized as follows. In Section II we 
present general formal expressions for the coherent and 
incoherent signals that allow to incorporate an arbitrary 
sequence of preparation pulses. The coherent signal is 
given by the time dependent polarization whereas the in- 
coherent is expressed in terms of transition amplitude. 
In Section III we skip the preparation stage (the shaded 
area in Fig. 2). We assume that the system has been 
prepared in a superposition state and calculate the spon- 
taneous light emission. In Section IV we apply the results 
of Section III to compute the spectra of cysteine induced 
by a stimulated X-ray Raman process. We further show 
that when the preparation is included as in Section II 
the coherent and incoherent signals represent difference 
frequency generation and fluorescence, respectively. 



II. COHERENT VS INCOHERENT 
SPONTANEOUS LIGHT SCATTERING SIGNALS 

Time-and-frequency gated spontaneous light emission 
signals is given by an integrated intensity of the electric 



field of an emitted photon. It may be further represented 
by the spectral and temporal overlap of a bare signal Wb 
and a detector spectrogram Wd [35] 



S(i,u>) = 



dt'^Wnii^t' ,u')W B {t' (1) 
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The detector offers a window of observation centered at 
time t and frequency Co. Even though i and Co can be 
varied independently the actual temporal <jt and spectral 
a w resolutions are not independent and satisfy cttc w > 1. 
This is guaranteed by Eq. (1). We shall assume a point- 
size detector, therefore omitting all effects of spatial res- 
olution only retaining temporal and spectral gating [35]. 
Eq. (1) was introduced originally as signal processing 
for the classical field [36] but it more broadly applies to 
quantum fields [35]. The bare signal Wb contains all 
of the relevant information about the externally driven 
molecules and is given by a correlation function of mat- 
ter that further includes a sum over the detected modes. 
In Appendix A we present superoperator expressions for 
W B of a system prepared in an arbitrary non stationary 
state without treating the preparation explicitly. The 
bare signal (A6) represented by the closed path time- 
loop diagram shown in Fig. 2a, b is given by time-ordered 
Green's functions for superoperators (TIGERS) [32]. It 
contains a double summation over molecular positions r a 
and rp and contains two types of contributions: incoher- 
ent emission terms with a = f3 (Fig. 2a) and coherent 
emission (Fig. 2b). 

We start with the coherent emission which is induced 
by pairs of molecules labeled a and j3. Here we can factor- 
ize the matter correlation in (A6) (TVj* (t{)V^ (h)) -> 

(Vl (*i))<yR t (*i)>- w b is then factorized into a product 
of amplitudes (see Appendices B, C) 

N N ~ l ikn-CR^-Rc) 



li JU a=l /3=1 a p 

x p( a -^(t' - t/2 - R a /c)P^^*{t' + t/2 - Rp/c). 

(2) 

Here the polarization of molecule a in the cartesian di- 
rection fi is given by 

$>aa ((&6|7%x(t)cxp (~ f H'_{T)dT) \aa)) a , 
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where p aa represents the initial equilibrium density ma- 
trix, T is responsible fo time ordering of superoperators. 
When using a purely temporal gate, Ff{uj,uj) = 1 and 
we get W D {Cj,t;t,T) = S(T)F*(t + t/2, - r/2,t), 
where W D {Q,t;t, T ) = ^e iuT W D (Q,i;t,uj). For an 
ideal time gate \F t {t,t)\ 2 =5(t-t) and Eqs. (1) and (2) 
result in the time-resolved signal 



S(t) = \E(t)\'< 
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where A is a normalization constant and the double- 
dot represents the second time derivative. In the op- 
posite, pure frequency, gating, i.e. F t (t,t) — 1, and 
the frequency gate is very narrow, such that Ff(t,cj) = 



v_y r —iujt- 



T*0(i) at 7 -> 0, then W D (u>, t; w, r) = e 



We then obtain the frequency-resolved signal 

S(Q) = \E(w)\ 2 (6) 
where E(lo) = dtE(t)e wt . Using Eq. (2) we obtain 
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It follows from the Larmor formula Eqs. (5) or (7) 
that only when the measurement is either solely time 
or frequency-gated, can the coherent signal be expressed 
in terms of the modulus square of a matter amplitude. 
The more general time and frequency gated signals (Al), 
(A6) and (2) do not have this form. 

We next turn to the incoherent signal (Fig. 2a) orig- 
inating from processes in which both signal-mode inter- 
actions occur with the same molecule a = j3. This sig- 
nal is not given by a modulus square of the polarization 
amplitude. However by expanding it in the molecular 
eigenstates we can express it in terms of other type of 
amplitudes. After a bit of algebra (see Appendix C for 
details) we obtain 
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represents a matter transition amplitude of molecule a 
starting with initial state a and reaching the final state 
b at time t. 

For a pure time domain gating (4) , Eq. (8) gives 
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Pure frequency gating (6) yields 
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Eqs. (8) - (11) indicate that the incoherent signal may 
not be generally recast in the form of square of the total 
polarization alone Eq. (3). Rather, it depends on a dif- 
ferent set of matter quantities (9), transition amplitudes, 
that represent the possible quantum pathways of matter 
from state a to state b. This calculation does not require 
superoperators since the amplitudes represent the evo- 
lution of the bra or ket separately and can be recast in 
Hilbert space. The polarization (3) on the other hand is 
governed by evolution of both bra and ket and therefore 
can only be defined using superoperators. 

The absence of spatial phase factors in Eq. (8) in- 
dicates that the incoherent emission is isotropic and is 
independent on the wave vector of the incoming pulses. 
The single summation over molecules yields linear ~ A^- 
scaling with the number of active molecules. The coher- 
ent signal (2) in contrast scales as N(N — 1). This sig- 
nal is directional and propagates according to the phase 
matching condition k s = k„. Clearly, the coherent and 
incoherent signals can be always distinguished by their 
N(N — 1) vs iV-scaling and directionality vs isotropic 
emission. A different type of size scaling exists in hyper- 
Rayleigh techniques. Consider a sample of TV particles 
each composed of M molecules. Assuming coherent in- 
traparticle and incoherent interparticle scattering the sig- 
nal will then scale as NM(M — 1) and will appear to be 
linearly scaling in the number of particles. iV-scaling can 
be similarly obtained if each molecule a is near some set 
of N' neighbor molecules /3 such that |k s | • |r Q — r^| <C 1, 
then it will emit in phase with this set of molecules, thus 
lending an overall factor of NN' . Incoherent sum fre- 
quency generation is an example of hyper-Rayleigh sig- 
nal [39-41]. This signal scales linearly with the number 
of molecules similar to the incoherent fluorescence. 

Classically, spontaneous light emission from a non- 
stationary state is calculated by the Larmor formula (see 
Appendix D). All the necessary information is then con- 
tained in the charge density of the accelerating parti- 
cle and the radiation power is related to the second or- 
der time derivative of the macroscopic polarization. The 
classical derivation apparently applied even to a single 
molecule and the result is then simply multiplied by the 
number of molecules. This yields linear TV-scaling. It fol- 
lows from our analysis that the classical formula is valid 
only for the coherent signal from a mesoscopic volume 
and may not be applied to a single molecule. Light scat- 
tering from a single molecule is incoherent and provides 
new information about quantum pathways of matter. It 
may not be recast using the charge density alone. 

To further evaluate the sum over molecules in Eqs. 
(2) - (11) we first consider small samples (compared 
to the optical wavelength) when retardation effects 
can be neglected. In this case both coherent and 
incoherent signals are independent of the direction of 
the incoming pulses. Straightforward analysis yields 
for the incoherent signal ^2 a Ra 2 ~ A^/i? 2 , where R c 
is characteristic size of the sample. Similarly since 
|k„| • |r a — r^l <C 1 the coherent component yields 



EaE^a^r 1 ~ NN'/R 2 C . Thus, both signals scale 
linearly with molecular density and inverse proportional 
to the square of the sample size. For extended samples, 
when retardation becomes important, the coherent 
signal yields a directional emission with iV 2 -scalmg: 

J2 e i(uiR a /c-k n -R a )-i(u'Rff/c-k n -Rp)^ a j^^-l ^ 

N(N - i)5(u/c - k„)5(u'/c - k n ). 

When retardation is neglected, the coherent signal (5) 
can be recast in the Larmor form. The most striking 
limitation is that this form may not be applied to the 
single molecule shown in Eq. (10); it requires 

pairs of molecules. The conventional classical derivations 
define the signal for a single molecule and then multiply 
it by the number of molecules. We showed that the N- 
scaling can be realized in the coarse grained system with 
mesoscopic grains via hyper-Rayleigh scattering where 
the short-range coherence makes the signal to scale N 
and not iV 2 . 



III. LIGHT SCATTERING OFF A 
NON-STATIONARY STATE 

We now calculate the incoherent and coherent signals 
in a multilevel system prepared in an arbitrary super- 
position state. In Section II, the evolution governed by 
H'_(t) and H'(t) in Eqs. (3) and (9), respectively caused 
by an arbitrary sequence of pulses which occurs prior 
to the last emission constitutes the preparation process 
that leaves the system in the superposition state. Here 
we ignore the details of the preparation (the shaded area 
in Fig. 2) and start our analysis after that time period 
where the system has been prepared in non stationary 
superposition state: p CC ' in Fig. 2a and two molecules in 
p cb and pbc'i respectively, in Fig. 2b. 



A. The coherent signal 

Starting with Eq. (3) and omitting the superscript 
a. For clarity we can recast the v-th component of the 
polarization in the form 

pM{t)=Tr[V v p(t)]. (12) 

After multiple interactions with incoming pulses the mat- 
ter density matrix Fig. 2b is p cb . We further define 
p^ = {b\V„\c}. We can then write: 



P M (t) = Yp b J c U b ) e-^ t . 



b,c 



(13) 



Assuming that molecules a, j3 are identical and neglect- 
ing retardation we obtain the following compact form for 
the bare spectrogram: 

W B ,coh(t',U>') 

= N(N ~ 1)A 2 / dTe- iu ' T P{t' ~ r/2)P*(i' + t/2), 
Jo 

(14) 
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FIG. 3. Various valence level-schemes considered in this pa- 
per, (a) - a general level scheme with arbitrary transition 
dipoles. (6) - a two band model with only interband transi- 
tion dipole, (c) - two bands with a single excited state, (d) - 
two bands with a single ground state. 



where we have further averaged over the random dipole 
orientation £ w P^ (t^P^* (t 2 ) ~ ipfa^fa) yield- 
ing A = A/y/3. Since P(t) only depends on a single 
interaction with the signal mode, the coherent signal gov- 
erned by the expectation value of the signal field is finite 
and the molecule returns to state \b)(b\ at the end of the 
process. In the case of ideal time-gating, we recover the 
Larmor formula: 



S coh (i)=N(N-l)^\P(t)\\ 



where 



P(t) = $>6c 



J bcV cb 



(15) 



(16) 



For purely frequency gating, we similarly get: 

Seoh(u) =NA 2 J2^LWb\ 2 \Pbc\ 2 S(cj~uj bc ). (17) 

b,c 



B. The incoherent signal 

To calculate the a — f3 terms in Eq.(A6) (see Fig. 
2a) we must evaluate the following matter correlation 
function: 
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where the state of the system prior to the emission is 
given by the nonstationary density matrix element p c > c 
and the emission brings it to the population state \b)(b\. 
Combining this with Eqs. (8) gives: 



W B ,inc(t' ',£*/) = nA 2 N ^ ^be^be'PcbP* d bPe 



b,c,c' 
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The incoherent bare spectrogram (19) may not be recast 
in the Larmor form (Eq. (14)). The time gated incoher- 
ent signal is 
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E ^bc^bc'P-cbPl'bPc'ce 1 ^ 
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(20) 



The frequency-gated signal is: 

S mc (u>) = NA 2 Y uj^ 2 cbPcc S(uj - uibc) (21) 



b,c 



In contrast with the coherent signal, we see that an 
excited-state population rather than a coherence in the 
relevant transitions is required to produce a signal. This 
is because all interactions are with a single molecule and 
the radiation signal mode (initially in a vacuum state) 
must be brought to a population one photon state to 
generate a signal. 

We now apply our results a level scheme composed of 
two manifolds of states with interband dipole elements 
as shown in Fig. 3b. Eqs. (14) - (19) which gives: 



W BOiCoh (t',uj') =irNA 2 
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(22) 



WWc(tV) = irNA 2 Y,^lgPgePge>Pee'e^-^^6 U , ( 23 ) 

gee' 



where the subscript "0" signifies the two-band model. 
In the case of time-and-frequency gating, various tran- 
sitions between ground and singly-excited states become 
coupled. A time-gated measurement (4) yields 



S ,eoh{t) = N- 



E 



"egVgePgee 



(24) 



S ,inc(t) = N—^2^2ojl g ojl g fj, ge p e > e e> {oj °'9 
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(25) 



For an ideal frequency gating (6), the signal (22) - (23) 
reads 

S 0tCoh {w) = NA 2 Yw%\fi ge \ 2 \p ge \ 2 S(w - w eg ) (26) 
eg 



S ,inc(u) = NA 2 ^2ujj g \ Pge \ 2 p ee S(uj - uj eg ) (27) 



Eqs. (26) and (27) are very similar, both given in 
the form of the Fermi golden rule for spectrally well sep- 
arated states. However, this is no longer the case for 
time-domain measurements. In Eq. (24), the initial and 
final states enter in the same way (the summations are 
fully interchangeable and all transitions add at the am- 
plitude level) while in Eq. (25), a trace is taken over final 
states after adding the amplitudes for the transitions to a 
given final state. The consequences of this difference are 
most readily seen by comparison of the two level schemes 
shown in Fig. 3. In case (c), there is a single excited state 
and a manifold of ground states. In case (d), the situa- 
tion is reversed and there is a manifold of excited states 
and a single ground state. Both cases will generate a 
beating term in the coherent signal. In contrast, the in- 
coherent signal (25) gives a beating term for case (d) but 
not for case (c). V and A type three-level model of atoms 
are commonly discussed in quantum optics [29, 42]; these 
correspond to our cases (d) and (c) respectively. The to- 
tal signal is S = S co h + S inc . S co h only has beats for (d). 
The semiclassical approach only gives S co h- 
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FIG. 4. (Color online) The stimulated X-ray Raman induced 
NLS process. Straight arrows correspond to interactions with 
X-ray pulses, wavy arrows represent spontaneous emission, 
(o) loop diagram for the bare SRIF signal (b) the bare co- 
herent DFG signal, (c) The level scheme used for cysteine is 
composed of a ground state go and 50 valence excited states 
gi with energies between ~5.75eV and ~11.5eV, / are core 
excited states. 



IV. DIFFERENCE FREQUENCY GENERATION 
VS FLUORESCENCE INDUCED BY A 
STIMULATED X-RAY RAMAN PROCESS 



We have simulated these signals for cysteine - an amino 
acid which serves as an important structural unit in con- 
necting different regions of proteins by disulfide bonds. 
Following excitation by stimulated X-ray Raman process 
resonantly tuned to either the nitrogen, oxygen or the 
sulfur K-edges. The signal is given by the diagrams of 
Fig. 4 and involves six radiation/matter interactions. In 
this case the incoherent signal is SRfF the coherent sig- 
nal is DFG as depicted in Fig. 4c. An attosecond X-ray 
pulse excites the core transition to state / by absorbing 
a ki photon which is then de-excited with — k2 photon, 
leaving the system in a superposition of valence states 
gi and 32 with amplitudes given by Eq. (A7) of Ref. 
[24]. This results in a weakly-excited, pure-state g%. We 
thus expect the incoherent and coherent signals to coin- 
cide for the 2 O go transition. The relevant diagrams 
are shown in Fig. 4a, b. For a two state model system 
prepared in a pure-state the state vector reads 



IV'o) = K g \g) +^K e |e), 



(28) 



where expressions for Kj, j = e,g in Eq. (28) are given in 
Ref. [24] . In the limit of weak excitation most of the pop- 
ulation is in the ground state \k 9 \ 2 ~ 1 and \K e \ 2 <C lVe 



(the excitation is perturbative) . Thus, Eq. (17) reads: 
S coh (Cu) = NA 2 J2"lWb\ 2 \Kb\ 2 K\ 2 5(u-u; bc ) (29) 



Since |K e | 2 <C 1 and 
yields 



1 the leading contribution 



S coh {u) =NA 2 J2"$M 2 \Kc\ 2 S{u>-Uge), (30) 

c 

where we assumed b — g. This is clearly a subset of the 
incoherent signal (Eq. (21)) which, under the assumption 
of a pure state, is given by: 

S inc (Q) = NA 2 J2^bcWb\ 2 \K b \ 2 S(Q-u bc ) (31) 

b,c 

Thus, for a weakly excited system in a pure state, the 
incoherent and coherent signals coincide for transitions 
from the excited states to the ground state. When the 
excited-state manifold bandwidth is smaller than the 
band gap, these transitions appear in the high-energy 
part of the emission spectra (see Section V). The intra- 
band transitions within the excited state manifold then 
appear in the red part of the spectra and only show up 
in the incoherent signal. 

The optimized geometry was obtained with the quan- 
tum chemistry package Gaussian09 [43] at the B3LYP[33, 
44]/6-311G** level. Core excited states were calculated 
with restricted excitation window time-dependent den- 
sity functional theory (REW-TDDFT) [34] implemented 
in NWChem code [45]. Core and valence energy lev- 
els and transition dipoles moments were calculated with 
a locally modified version of NWChem code at the 
CAM-B3LYP[46]/6-311G** level of theory and within 
the Tamm-Dancoff approximation [47]. Additional com- 
putational details are given in Ref. [25]. The calculated 
frequency gated SRIF (21) and DFG (17) signals ob- 
tained by excitation of the nitrogen, oxygen and sulfur K 
edges are shown in the top and bottom rows respectively 
of Fig. 5. 

We first note that the SRIF and DFG signals coincide 
in the high-frequency regime. The lowest excited state 
has an energy ~ 5.74 eV while the excited state manifold 
bandwidth is ~ 5.70 eV. These high-energy transitions 
are from the various excited states to the ground state. 
The low-energy us < 5 eV features that only appear in 
the incoherent spectra represent transitions between ex- 
cited states gi,2 f-> 33. The frequency gated coherent 
and incoherent signals are very different. However the 
corresponding time gated signals shown in Fig. 6 are less 
distinct. 

Fig. 7 depicts the time-and-frequency gated spectro- 
grams of the SRIF signal computed by Eq. (14) of Ref. 
[35] with gaussian gating functions 
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FIG. 5. (Color online) SRIF (Top) and DFG (Bottom) NLS signals of cysteine following stimulated Raman excitation by 
an X-ray pulse of width 14.2 eV and central frequency tuned to the nitrogen (~ 404.4 eV), oxygen (~ 532.2 eV), or sulfur 
(~ 2473.5 eV) K-edges, as indicated. The two signals coincide for gi : 2 go transitions. Transitions between excited states 
gi,2 o g3, uj < 5 eV do not show up in the coherent signal. F ~0.04 eV is used for the gating bandwidth. 
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FIG. 6. (Color online) Time resolved SRIF - solid red line 
and DFG - dotted blue line signals are compared for X-ray 
pulses resonant with S, O and N in cysteine. F ~0.04 eV is 
used for the gating bandwidth. 



where <r w and ot are the corresponding bandwidths of 
the frequency and time gates, respectively. These obey 
Fourier uncertainty relation a^UT > 1- Some transitions 



in Fig. 7a do not vary with time while others beat. Beat- 
ing occurs because of interference with nearby transitions 
and thus indicates a higher density of transitions. The 
beating period gives the interval between the interfering 
states a, b via \u> a b\ = T ■ The gating parameters de- 
termine the maximum energy difference (and thus the 
minimum beating period) between interfering states via 

A = a L0 + — , (33) 

where A is the maximum spacing between interfer- 
ing states. Note that these resolutions always satisfy 
&u&t > 1 [36]. One can thus begin by applying a 
nearly-pure frequency-gating to eliminate interferences 
and then gradually widen the gating parameters in the 
frequency-domain to reveal the density of transitions at 
various parts of the spectra. A clear demonstration is the 
prominently-beating transitions with period ~ 120 fs at 
~ 9 eV which is visible (to varying degrees) in all three 
2D figures. An optimal balance for observing this inter- 
ference is found in Fig. 7a in which the beating period 
(and thus the state seperation) can be well-approximated 
by visual inspection. In Fig. 7b A is so large that the 
interference of interest is masked by the myriad of other 
possible interferences while in Fig. 7c, A is small enough 
to almost entirely suppress the ~ 9 eV beating. The 
coherent DFG signal coincides with SRIF in the regime 
Q > 5 eV and vanishes otherwise. It is thus, not plotted 
separately. 

We note, that the choice T — 0.04 eV in Figs. 5, 6 
and 8 represents an experimental resolution. The life- 
time of the valence states, are of the order of ns and will 
give broadening r ~ [i eV. The detection parameters are 
chosen to capture the specific beating signal between var- 
ious transitions that occur on much shorter time scale of 




w (eV) 

FIG. 7. (Color online) Time and frequency gated SRIF sig- 
nals of cysteine (Eq. (1)). The gating parameters in atomic 
units a T = 1000, cr u = 0.001 - (a), a T = 100, cr^ = 0.01 - (6), 
and or = 2000, — 0.005 - (c). These permit all transitions 
within ~ 0.054 eV, 0.54eV and 0.1 eV of each other to inter- 
fere, respectively. Top panel has optimal gating parameters 
and reveals both areas of low transition density (where the 
intensities are time independent) and areas of higher densi- 
ties (where beats develop as a result of interference between 
transitions). Note the particularly prominent beating near 9 
eV with an approximate period of 120 fs. The actual distance 
between these states is ~ 0.0126 Ha well within the allowed 
detection bandwidth. Middle panel has low frequency resolu- 
tion but rather high temporal resolution that results in clear 
beating signal. The bottom panel has high spectral but low 
temporal resolution, that results in the suppression of for- 
merly prominent beating at ~ 9 eV. Further narrowing the 
frequency-domain gating widths eliminates it altogether. 



tens of fs. The low quantum yield at short times poses a 
challenge for the detection. 



V. DISCUSSION 

Thanks to the ~ iV 2 versus ~ TV-scaling it is not hard 
to separate the coherent signal since it is directed and 
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FIG. 8. (Color online) Frequency-resolved SRIF - green (light 
gray) line and DFG - red (dark gray) line following a stimu- 
lated Raman excitation by two pulses of varying width. In all 
three panels, the first pulse is tuned to the transition from the 
bottom of the nitrogen core band to the ground state (~ 388 
eV) and has a width of 1.36 eV. The second pulse is tuned 
to the transition from the bottom of the nitrogen core band 
to the top of the valence band (~ 374 eV). The second pulse 
is gradually broadened to allow for more valence excitations 
prior to the NLS ~1.36 eV (bottom), ~3eV (middle), and ~6 
eV (top). F ~0.04 eV is used for the gating bandwidth. 



generally dominates the incoherent and hyper-Rayleigh 
coherent signals. However, it is not generally as simple to 
distinguish between the incoherent and hyper-Rayleigh 
contribution to the coherent signal since both are roughly 
isotropic and ~ ./V-scaling. In our model of cysteine, the 
valence manifold bandwidth is slightly larger than the 
lowest valence excitation and so the incoherent and co- 
herent signals can be easily distinguished in the frequency 
domain, as shown in Section IV. 

When the valence state manifold bandwidth is compa- 
rable to the band gap, both inter- and intra-band transi- 
tions overlap spectrally and one needs to identify another 
procedure for separating the SRIF and DFG contribu- 
tions. A single X-ray pulse cannot accomplish this goal 
since any modification of its bandwidth will affect both 
the absorption and the emission profile. However, this 
can be achieved by using two X-ray pulses with wave 
vectors k! and — k 2 for the system preparation via stim- 
ulated Raman excitation. Both pulses must have narrow 
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bandwidth compared to the g manifold in order to excite 
states selectively. The first pulse is tuned to the bot- 
tom of the core-excitation manifold. The second pulse 
is tuned to the transition from the bottom of the core- 
excitation manifold to the top of the valence band. This 
combination of pulse parameters can selectively excite 
populations and coherences of the high-energy valence 
states. Since the excitation is weak, the square of the co- 
herences ground and excited states are comparable to the 
valence populations but the coherences between valence 
states are lower order. Thus, the coherences between 
valence states created by this two-pulse process do not 
contribute significantly to the signal. The overall signal 
will then be dominated by the incoherent contribution, in 
particular the transitions from high energy valence states 
to lower energy valence states and the ground state. In- 
creasing the bandwidth of pulse — k 2 will result in the 
excitation of larger number of valence states which gives 
rise to a small coherent contribution as well as more tran- 
sitions to appear in the incoherent spectra. Finally, as 
the second pulse is broadened so that it excites the entire 
valence band, the incoherent and coherent signals again 
overlap in the regime of transitions to the ground state. 
This behavior is illustrated in Fig. 8. 

In summary, by working in the joint field-matter space 
we developed QED expressions for the NLS signal and 
show that may not be solely expressed in terms of the 
macroscopic polarization as suggested by the semiclassi- 
cal approach. We find an isotropic hyper-Rayleigh and 
directional coherent component that can be described by 
the Larmor formula and an incoherent component that 
governs the fluorescence. The microscopic calculation re- 
veals interferences of the quantum matter pathways. 
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Appendix A: Microscopic Calculation of 
Spontaneous Light Emission 

Simultaneously time-and-frequency gated spontaneous 
light emission signals may be described by the spectral 
and temporal overlap of a bare signal and a detector spec- 
trogram (Eq. (1)). The detector spectrogram Wd is an 
ordinary function of the gating time and frequency pa- 
rameters which are characterized by standard deviations 
of the time and frequency gating gt and <r w , respectively. 
The structure of Wd guarantees that these always satisfy 



the Fourier uncertainty (j^/jt > 1. The bare signal con- 
tains all of the relevant information about the molecules. 

In order to maintain the bookkeeping of all interac- 
tions and develop a perturbative expansion for signals we 
adopt superoperator notation. With each ordinary oper- 
ator O we associate a pair of superoperators [32] "left" 
6 L X = OX, "right" R X = XO, and the combination 
0_ = Ol — Or. The bare spectrogram Wb in the gated 
photon counting signal (1) is given in terms of superop- 
erators as: 

/•OO 

W B {t',J) = / dre- tul ' T x 
Jo 

{TE R (t' + r/2, r D )± L (t' - r/2, v D )e~^ /~~ A '- {T)dT ). 

(Al) 

The positive frequency part of the electric field operator 
is given by 




(A2) 

where e"*)(k) is the unit electric polarization vector of 
mode (k s ,/i), fj, being the index of polarization, uj s = 
c|k s |, c is speed of light, Q is quantization volume. The 
polarization degrees of freedom are necessary to restrict 
the density of radiation modes and recover the A~ 4 de- 
pendence typical for Rayleigh scattering. The Hamilto- 
nian superoperator in the interaction picture under the 
rotating-wave approximation (RWA) is given by: 

H' q (t) = J dr¥,\{t,v)\ q {t,v) + H.c, q = L,R, (A3) 

where V(t, r) = J2 a V Q (t)<5(r — r a ) is a matter operator 
representing the lowering (exciton annihilation) part of 
the dipole coupling and a runs over molecules in the sam- 
ple located at r a , rp. The operator T maintains positive 
time ordering of superoperators, and is a key bookkeep- 
ing device. It is defined as follows: 

t.e,(*i)4'(*2) = e(h - t 2 )4(*i)4'(*2) 

+ e(t 2 -t 1 )E q ,{t 2 )E q {t 1 ), (A4) 

where 6(t) is the Hcavisidc step function. Note that the 
interaction Hamiltonian (A3) corresponds to the interac- 
tion of the matter with the spontaneously emitted pho- 
ton between two electronic states. This energy is in VUV 
range and is not related to the preparation of the system 
that involves an X-ray Raman process with high energy 
photons. Therefore, the dipole approximation for field- 
matter interaction is justified. 

We shall calculate the time dependent bare signal (Al) 
for a collection of molecules a, [3 in the interaction pic- 
ture. We consider a signal governed by the spatial phase 
factor e lk "' r , where k„ is a combination of incoming wave 
vectors characteristic to the desired signal. This phase 
factor represents a paraxial approximation, such that 
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multiple frequency components propagate with the same 
wave vector. The bare signal represented by the closed 
path time-loop diagram shown in Fig. 2a,b. We first 
expand the time ordered exponent in Eq.(Al) to second 



order - one interaction with the bra and one with the kct 
- and factorize the detected field and matter correlation 
functions 



1 N r t'+r/2 rt'-r/2 

/ dtx ' ,W rT\'<< ,,i 11- *1 ,\/T i;< , ,' ... .. /•) ,. -;r- -../■) ... ,r ( " 'hi' .. \ r" : 



h 2 



a,0 vjs' 



Wx(TV? L &)V$\{tMTE™{i! + T/2,r D )E^(t>-T/2,r D )Ell>\t[,r a )ElV R (t ll re)}. 

(A5) 



Since Eh'l is initially in the vacuum state, the field correlation function factorizes using (A2) and the bare signal (Al) 
is given by 



W B {t',Lj' 



r oo r t'+r/2 r t'-r/2 

J — oo J— oo 



k 3: k 



N 



x ^ el (k s - k „).R,-,( k3 ,- k „).R Q £ e M(k s )eW(k s )e^')(k s 0^^^ (A6) 

a,/3 v-^v-v' 

where R a = r Q — and yd is the position of the detector. Eq. (A6) contains explicitly multiple pairs of radiation 
modes k s and k s / and acts in the joint field plus matter space. It takes into account all field matter interactions that 
lead to the emission of the detected field modes. 



Appendix B: Superoperator representation of coherent and incoherent signals 

For the single-molecule (incoherent) signal r a = rg and thus, the signal reads 



W Binc (t',J) 



k a ,k 



N 



, Jo 



t'+r/2 ft'-r/2 



]T / dre-^' T J2^ K ' yRa Mi dt\u s J s e^ {t ' + h^)-^'-h2-h) 



£ e(^(k s )eW(k s )e^)(k s 0eW(k s 0(TK i ^)^ iJ (t 1 )e- 



-iS-ooH'-(T)dT 



(Bl) 



In the absence of dissipation (unitary evolution) we can further factorize the matter correlation function as 
{TVMV^ih)) = ^ Paa ((aa\V^ R (t 1 )\ab)){(ba\V vL (t' 1 )\aa)) = 



a.b 



J2 Paa{W + T/2)a\V uL (t' 1 )T+ exp 

a.b 



— J^ H' L (T)dT \aa))((aa\vS R (h)T-exp H> R (T)dT 



W+T/2)a)), 
(B2) 



where we denote ((ba\0\aa)) = Ti[\a) (b\0\a) (a\] and b is the final state of the system. We assume that system is 
initially in the pure state a described by equilibrium density matrix p aa . We next define the transition amplitude 



ioj s (t— t^) — iuj a bt— i(\t s — k{^| )-R Q 



x {{b{t)a\V vL {t' 1 )Te-x. V i- 1 - J ' H' L (T)dTj \aa)) 



(B3) 
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Since all interactions are from the left (L), we can also write the transition amplitude using ordinary operators in 
Hilbert space 



^r ) (*)=-E eH ( k ^(k s ) 2 ^/ i 

x (6(t)|K(*i)Texp f - l - p H'(T)dTj \a) a (B4) 



-iuj s (t— — ioj a bt— «(k s — k{^| )-R Q 

— i - 

f*'i 



This gives for the bare signal (Bl) 



o,6,M a=1 

x ^ r ) ( ^_ r/2) f(^)*( t ' +T /2). (B5) 

In the limit of pure time-resolved measurement (4) signal 
(B5) transforms into 

JV 

W Blnc (t) Paalf^Ht)] 2 , (B6) 

ci.b.fi a = l 

Similarly the pure frequency-resolved signal (6) signal 
(B5) yields 

JV 

W Binc (u>) Paa^m 2 , (B7) 

a,b,fia — 1 I 



where f£'»\w) = dtf^ (t)e iut . The spatial 

phase factors in Eq. (A6) indicates that the incoher- 
ent emission occurs in all directions and is independent 
of the wave vector of the incoming pulses as expected. 
A single summation over the molecule yields the linear 
scaling with respect to the number of molecules. 

We now turn to coherent emission contribution of a =/= 
(3. Since interactions with different molecules arc not 
time-ordered, we may factorize the matter correlation in 
(A6) (TVZWVlHh)) -+ {Vm))(ViHh)). Thus, the 
two-molecule (coherent) signal can be separated into a 
long-range and a short-range component (as described in 
[48]). Using Eqs. (B2) - (B3) the coherent part of the 
signal (A6) reads 

W Bcoh {t',J) = 

JV JV-1 oo 

E E E / dTe- iw ' T P^ (t 1 - t/2)P {!3 ^ {t 1 + r/2) 

fj, a =l /3=1 ^° 

(B8) 

where 



P< a -">(t) - -iEe W (k s ) £ W(k 8 )^ f dt ' ie -<«-.(*-*i)-*-.»t-i(k.-k {B} ).R 
x Ep aa ((66|K L (i / 1 )Texp (^- % - J 1 H' L {T)dT^j \aa)) a . 



Here p aa represents the initial equilibrium density ma- 
trix. Note that in contrast with Eq. (B3) which repre- 
sents the transition amplitude and thus can be recast in 
Hilbert space, (B9) is related to polarization and there- 
fore can be written in Liouville space only. In the limit 
of pure time-resolved measurement (4) signal (B8) trans- 
forms into 

JV JV-1 

w Bcoh {t) = E E E p (a ^(t)p {M (t), (bio) 

H a=l 0=1 

Similarly the pure frequency-resolved signal (6) signal 



(B9) 



I 

(B8) yields 

N JV-1 

w Bcoh {u>) - E E E p m (")p {M ("), (bii) 

H a=l /3=1 



Appendix C: Microscopic derivation of the Larmor 
formula 

In order to calculate the signal (A6) we start with am- 
plitude expression (B3) . We first evaluate the summation 
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over the modes. In the continuum limit 

Recalling that k s and polarization vectors e^(k s ), v = 
1, 2 form a set of mutually perpendicular unit vectors it 
follows that 



Taking into account (B5) we therefore obtain (8). Simi- 
larly we derive for coherent signal 

P^\t) = -^ e ^£ld 2 t P^ a \t- R a /c). (C8) 

Thus, coherent signal Eq. (B8) can be recast in the form 
of (2). 

Appendix D: Semiclassical theory of emission 
detection from a radiating dipole 



e M(k s )eM(k s ) = S^-kX- 



(C2) 



We further evaluate the momentum integral using dk s 
k 2 dk s dQ. The angular integration yields 



47T 



dft(5 M „ - k^k v )e 



ik s R 



= p(-V 2 V + V,-V,) Sin( |^. (C3) 



Eq. (B3) then yields 



R 



fir\t) = e ^Y,it T ^ ) ^ ( 

V 

sin(co s R a /c) 1 



/ 



x / duo s {-V z 8^ + V M • V„)- 



i? a w — cj s + it ' 
(C4) 



where we introduced a matter transition amplitude 
Ti, a {uj) that connects initial state a with final state b: 
^'"V) - n^-T^-^) and T b a f\t) is given by 
Eq. (9). We then note that 



u; s — lo — ie 



(C5) 



Taking into account that [37] 



(-v 2 v + v M - V„)- 



R 



jp[(6p, - 3RM(ikR - 1) + - R^R u )k 2 R 



2 D 2i e ikR 



(C6) 



and assuming the random orientation of the molecules: 
Rfj,R v = ^Sfiu we obtain 



According to classical electrodynamics, the electric 
field obeys the usual homogeneous wave equation derived 
from Maxwell equations: 

„ „ „, . 1 <9 2 E(r,t) 1 <9 2 P(r,i) 



e c 2 at 2 



Here P is the macroscopic polarization and e is the vac- 
uum permittivity. Since Eq. (Dl) is linear it applies also 
for a quantum field, where the electric field E and P be- 
come operators. Restricting ourselves to the transverse 
part in far field E = Ej_, this equation can be recast in 
frequency domain as 



2 2 

V 2 E(v,lo) - ^E{v,u) = ^P(r, W ) 
c e c 



(D2) 



The Green's function solution of Eq. (D2) in infinite 
space for a single-point dipole molecule at r = 



E(t,u>) 



1 P *fco|r| _ 

e c 2 47r|r| 



(D3) 



where k = uo jc. Latter can be recast in the time domain 

£;(r ' i) ^^4^^ ( °'*- |r|/c) (D4) 

Neglecting retardation effects the electric field entering 
the detector at position r G is given by 



E(r G ,t) = BP(t), 



(D5) 



where B = (47re c 2 |r G |)~ 1 and P(t) = P{0,t). We thus 
recover the Larmor formula for the time-resolved inten- 
sity of the signal 



S(t) = B 2 \P(t)\ 2 . 



(D6) 



Generally the signal will be affected by the detector ge- 
ometry, depending if an angle is used with a lens or if it 
enters in parallel. These details will affect the prefactor 
B and are not considered here. The bare signal can be 
finally calculated as 

poo 

W B {t',uj') = B 2 / dTe- lu> ' T P*{t' + T/2)P{t' -r/2). 
Jo 

(D7) 
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